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Disordered quantum gases under control 

Laurent Sanchez-Palencia 1 and Maciej Lewenstein 2 

When attempting to understand the role of disorder in condensed-matter physics, one faces severe experimental 
and theoretical difficulties and many questions are still open. Two of the most challenging ones, which have been 
debated for decades, concern the effect of disorder on superconductivity and quantum magnetism. Recent progress 
in ultracold atomic gases paves the way towards realization of versatile quantum simulators which will be useful 
to solve these questions. In addition, ultracold gases offer original situations and viewpoints, which open new 
perspectives to the field of disordered systems. 


P hase coherence and interference effects underlie many 
basic phenomena in mesoscopic physics, for instance 
electronic conduction 1 , magnetism 2 , superfluidity and 
superconductivity 3 , or the propagation of light and sound waves in 
inhomogeneous media 4 . Both also play central roles in high-precision 
devices such as interferometers, accelerometers and gyroscopes. In 
this respect, an important issue concerns the effects of disorder, i.e. of 
small random impurities, which cannot be completely avoided in real- 
life systems. A priori, one may expect that weak disorder slightly 
affects most physical systems and that averaging over the disorder 
smoothens possible effects. One may also expect that, in quantum 
systems, the spatial extension of wavefunctions leads to even weaker 
effects, via a kind of self-averaging. In fact, these naive ideas turn 
out to be wrong. Disorder often leads to subtle situations in which 
strong effects survive averaging over the disorder 4 , in particular in the 
quantum world. One of the most celebrated examples is Anderson 
localization 5 (AL). It is now understood that AL results from interfer¬ 
ence of the many paths associated to coherent multiple scattering from 
random impurities, yielding wavefunctions with exponentially decay¬ 
ing tails and absence of diffusion 6 . This strongly contrasts with the 
Drude-Boltzmann theory of classical transport, which predicts that in¬ 
coherent scattering induces diffusion 1 . 

Anderson localization was first introduced for non-interacting 
quantum particles to explain the absence of electronic conduction in 
certain dirty solids 5 , but remained elusive for matterwaves. It was re¬ 
alized later that it is actually ubiquitous in wave physics 4 , paving the 
way for the first observations of AL, using classical waves, e.g. light in 
diffusive media 7 - 8 and photonic crystals 9 - 10 , microwaves 11 and sound 
waves 12 . In condensed-matter physics, AL is now considered a fun¬ 
damental phenomenon underlying certain metal-insulator transitions, 
but complete theory of disordered solids should incorporate Coulomb 
interaction, the underlying crystal structure, interaction with phonons, 
and magnetic effects. Unfortunately, understanding the physics of even 
the simplest models including all ingredients poses severe difficulties 
and many issues are still unsolved or even controversial. The most 
challenging ones concern the interplay of disorder and inter-particle 
interactions, and spin-exchange couplings. 

Surprisingly enough, atomic physics offers new approaches to 
these issues. The field of ultracold atoms has been developing rapidly 
in the past decades, making it possible to produce, probe and ma¬ 
nipulate Bose 13,14 and Fermi 15 ' 16 gases with unprecedented versatil¬ 
ity, tunability and measurement possibilities (Box 1). Control in 
these systems is now such that ultracold atoms can realize quan¬ 
tum simulators 17,18 , i.e. platforms to investigate various fundamental 
models 19-22 . Landmark results have already been obtained, e.g. ob¬ 
servation of Mott insulators 23-25 , Tonks-Girardeau 26 ' 27 , Berezinskii- 
Kosterlitz-Thouless 28 physics, and magnetic-like exchange 29 - 30 . In¬ 
vestigation of Bose-Einstein condensates (BECs) in disordered 
potentials 31 - 32 has also emerged in a quest for direct signatures of AL 
of matter-waves. Joint theoretical 33-36 and experimental efforts 37 ^ 13 
made it possible and two groups succeeded recently in observing one¬ 


dimensional AL 44 - 45 . 

Prima facie, the discovery of this ‘Holy GraiT might mean the 
end of a quest. On the contrary, it is just a beginning as the two ex¬ 
periments of refs. 44, 45 open unprecedented paths to pursue many 
outstanding challenges in the field of disordered systems. Direct ex¬ 
tensions include studies of metal-insulator transitions in dimensions 
larger than one, and of the effect of weak interactions on localiza¬ 
tion, for which many questions are debated. For stronger interactions, 
single-particle localization is usually destroyed, but new concepts such 
as many-body Anderson localization 46-48 and Bose glass 49-52 provide 
original paradigms, which renew our understanding of these issues. 
Experiments on ultracold atoms with controlled disorder and con¬ 
trolled interactions can also be extended to other systems where dis¬ 
order plays important roles. For instance, combining spin exchange 
implementation 29,30 and disorder opens the route towards random 
field-induced order 53-55 and spin glasses 56-59 . These few examples il¬ 
lustrate all the promises of an emerging field, i.e. quantum gases in 
controlled disorder. In this paper, we review theoretical and exper¬ 
imental progress in this area and discuss perspectives that are now 
within our grasp. 

The nature of Anderson localization 

Localization, as introduced by P.W. Anderson in 1958, is strictly 
speaking a single-particle effect 5 . Consider the wavefunction of 
a free particle of mass m and energy E, in a d-dimensional quenched 
disordered potential L(r), which is solution of the Schrodinger equa¬ 
tion 

W- r) = --r—-<Kr) + VW(r). (1) 

2m 

While in free space, is an extended plane wave, it can be shown 
rigorously 60 - 61 that, in the presence of disorder, any solution with arbi¬ 
trary E is exponentially localized in ID, i.e. ln(|i/dL)|) ~ |z|/L loc , with 
localization length L loc (£) oc / B , where / B is the transport (Boltzmann) 
mean-free path. Eventhough Li oc often increases with E, it is striking 
that interference effects of multiply scattered waves is strong enough 
to profoundly affect even for very high energies. In 2D. the 
situation is similar 62 , but interference effects are weaker, and L loc oc 
/b exp(7rk/ B /2) where k = V2ra£/S would be the particle wavevector 
in free space. Hence L ]oc increases exponentially for k > 1// B , induc¬ 
ing a crossover from extended to localized states in finite-size systems. 
The situation differs dramatically in 3D where a proper phase tran¬ 
sition occurs at the so-called mobility edge k mob : While low-energy 
states with k < L mob are exponentially localized, those with k > k mob 
are extended. The exact features of the mobility edge are unknown, 
but approximately captured by the Ioffe-Regel criterion 63 , which ba¬ 
sically states that localization requires the coherence volume contain 
several scattering processes. In other words, coherence must survive 
on longer distances than the memory of the initial particle direction, 
thus yielding k mob ~ 1// B . 
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Direct observation of Anderson localization of matter- 
waves 

Observing AL of matter-waves requires several challenging condi¬ 
tions. First, one must use weak-enough disorder so that interference 
effects at the origin of AL dominate over classical trapping in po¬ 
tential minima. Second, one must eliminate all perturbations such 
as time-dependent fluctuations of the medium, or inter-particle inter¬ 
actions. Finally, one must demonstrate exponential localization, not 
only suppression of transport as it can also arise from classical trap¬ 
ping. While these conditions are very demanding in condensed-matter 
physics, they can be accurately fulfilled with ultracold atoms, using 
i) controlled disorder, ii) negligible interactions, iii) strong isolation 
from the environment, and iv) direct imaging of atomic density pro¬ 
files. This way, direct signatures of AL of non-interacting matter- 
waves were reported in refs. 44,45. As we shall see, these two ex¬ 
periments are complementary rather than similar because they signifi¬ 
cantly differ as regards both observation scheme and class of disorder. 

In ref. 44, a weakly interacting BEC is created in a trap, which is 
abruptly switched off at time t = 0. Then, the condensate expands in 
a guide and in the presence of disorder (Fig. 0a), created with opti¬ 
cal speckle (Box 2). This physics is captured by the Gross-Pitaevskii 
equation 

dilt Ti 2 V 2 , 

A— = ——ifr + V(r)i (/ + (2) 

which corresponds to Box 1 Hamiltonian (1) in meanfield regime. The 
dynamics of the BEC can be understood in a two-stage scheme 35,36 . 
First, it is dominated by interactions and the BEC expands, creat¬ 
ing a coherent wavefunction with a stationary momentum distribution, 
D(k) oc 1 - (fe£) 2 , where £ = fi/ -^4m^i is the initial healing length, 
which measures the initial interaction strength 13 . Second, once the 
expansion has strongly lowered the atomic density |ife(z)| 2 , the inter¬ 
action term vanishes and we are left with a superposition of (almost) 
non-interacting waves if/ k , the population of each is D{k). Then each 
i//k eventually localizes by interacting with the disordered potential, 
so that ln(|i/r t .(i;)|) ~ |z|/Li oc (fe), and the total BEC density reduces 
to 35,36 nBEc(z) — f dk D(fe)(|i//>(z)| 2 ). Direct imaging of the localized 
matter-wave reveals exponentially decaying tails 44 , with a localization 
length equal to that of a non-interacting particle with k = 1 /£ (Inset of 
Fig. |a). Hence, this experiment corresponds to a ‘transport scheme’, 
which probes AL of non-interacting particles with a wavevector con¬ 
trolled by the initial interaction, via 

In contrast, the experiment of ref. 45 uses to a ‘static scheme’. 
The interactions are switched off already in the trap via Feshbach res¬ 
onances, so that the gas is created in a superposition of a few (typically 
1 to 3) low-energy, single-particle eigenstates. They are subsequently 
imaged in situ, revealing exponentially decaying tails (Fig. 0b). It is 
worth noting that ref. 45 uses a ID quasi-periodic, incommensurate 
lattice (Box 2), thus realizing the celebrated Aubry-Andre model 64,65 

H = — ^ j(a)a, + h.c.) + ^ Acos(27r/?y + (f >) a^.cij (3) 

< 7,0 7 

i.e. Box 1 Eq. (2) with U = 0, Vj = A cos(2 nfij + (/>), and /3 an irra¬ 
tional number. Differently from the case of truly disordered potentials, 
there is a metal-insulator transition (mobility edge) in ID, which is 
theoretically expected at A/J = 2. 

These works open new horizons to further deepen our knowl¬ 
edge of AL in various directions. In ID, although all states are 
localized, subtle effects arise in correlated disorder, for instance in 
speckle potentials 35 . To lowest order in the disorder amplitude, Vr = 
sJ(V(z) 2 ), the Lyapunov exponent, y(fe) = 1 /Li oc (fe), can be calculated 
analytically 66 and one finds y(fe) oc (V(2k)V(-2k))/k 2 , enlightening 
the role of coherent second-order back-scattering, +k —» - k —» + k, 
in the localization process. Since the power spectrum of speckle po¬ 
tentials has a cut-off suc h that C2(2fe) = (V(2k)V(-2k)) = 0 for 
fe > fee (Box 2), one finds an abrupt change (effective mobility edge) 
in the fe-dependence of y for weak disorder 67,68 : While y(fe) ~ for 
fe < fee, higher-order scattering processes dominate for fe > fee and 
7(fe) ~ V*. 

In dimensions higher than one, the self-consistent theory of 
localization 69 allows one to calculate L loc and exhibits a mobility edge 



Figure 1 | Experimental observation of Anderson localization of mat- 
terwaves with Bose-Einstein condensates, a) Experiments of Institut 
d’Optique (coutesy of V. Josse and P. Bouyer): An interacting BEC expands 
in a tight ID guide (in red) in the presence of a speckle potential (in blue). 
The expansion stops in less than 500ms and the density profile of the con¬ 
densate is directly imaged (shown in orange-green; from the data of ref. 44). 
The column density, plotted in semi-logarithmic scale in the inset, shows a 
clear exponential decay characteristic of Anderson localization. The localiza¬ 
tion length L loc , extracted by fitting an exponential exp(-2[z|/L loc ) to the exper¬ 
imental profiles 44 , shows a good agreement with theoretical calculations 35,36 , 
b) Experiments of LENS (adapted from ref. 45 with permission of the au¬ 
thors): A non-interacting BEC is created in a combination of a harmonic trap 
and a 1D bichromatic lattice. The plot shows the exponent a of a fit of a 
function exp(-|(.r — x 0 )/l\ a ) to the tails of the condensate at equilibrium in the 
combined potential, versus the ratio of the disorder strength (A) to the site-to- 
site tunnelling rate ( J ). The onset of localization corresponds to the crossover 
to a -> 1 for A/J > 9. The inset shows a plot of the density profile of the 
condensate together with the fit for A/J = 15. 


in 3D. It is however known that it is not fully exact. Therefore, a ma¬ 
jor challenge for disordered, ultracold atoms is to extend the works of 
refs. 44,45 to two 70,71 and three 72 dimensions. Definitely, observing 
the 3D mobility edge would be a landmark result, which may stim¬ 
ulate further theoretical developments and drive new approaches by 
providing precise measurements of the mobility edge fe mo b and the cor¬ 
responding critical exponents, which are unknown. 

Interactions versus Anderson localization 

Another outstanding challenge is to understand how interactions affect 
localization, a question that has proved puzzling from the earliest times 
of AL 7 ’, and which is still debated. Common belief is that even weak 
interactions destoy localization. Different approaches however provide 
apparently contradicting answers in different transport schemes. For 
instance, recent numerical calculations 74,75 suggest that for expand¬ 
ing BECs, repulsive interactions destroy AL beyond a given thresh¬ 
old. Conversely, other recent results 76 predict that localization should 
persist even in the presence of interactions. Finally, in transmission 
experiments (which amount to throw a mono-kinetic wavepacket to a 
disordered region and measure transmission), perturbative calculations 
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Box 1 | Ultracold quantum gases. 


Creating and manipulating ultracold gases 

Ultracold quantum gases are dilute atomic systems that are cooled down to temperatures of the order of a few tens of nano-Kelvins and 
confined in immaterial traps using combinations of magnetic and optic fields 1416 . Owing to strong dilution, the prominent inter-particle 
interactions are two-body interactions while many-body interactions can often be ignored. At ultra-low temperatures, s-wave scattering 
dominates and the interaction is accurately described by a contact pseudo-potential 1315 . In the general case of mixtures of atoms in 
different species (or different internal states), the physics is thus governed by the Hamiltonian 


" = Z/ dr V ^ (r) 


ft 2 V 2 
2 m„ 


+ ^(r) 


'Kr(r) 




(1) 


where TV and m a are the field operator and the mass of an atom of species a. The first integral in Box 1 Eq. (1) represents the single¬ 
particle Hamiltonian where the potential W(r) is controlled by the configuration of the magnetic and/or optic fields (Box 1 Fig la). In 
most cases, it is nearly a harmonic trap 14,16 (VV(jc, y, z) = Zrep-.i^) "t (r <jj 2 ,/" 2 /2), the anisotropy of which can be adjusted in experiments. 
For instance, making it strongly anisotropic offers the possibility to produce one- 105,107 or two- 28 dimensional gases. Another useful 
possibility is to create a guide for the atoms using a strongly focused laser beam 108 . The second integral in Box 1 Eq. (1) represents 
the interaction operator where g, ur i is the coupling constant for interacting atoms of same or different species (go-,a' > 0 and g^ < 0 
correspond to repulsive and attractive interactions, respectively). Interestingly, the value and the sign of g CT(r - can be controlled in quantum 
gases using Feshbach resonances 22 . 


Optical lattices 

Considering different limits of Hamiltonian (1) allows one to design various models initially introduced in the context of condensed-matter 
physics, but here in a controlled way. One important example is that of optical lattices, which are produced from the interference pattern 
of several laser beams 20-22 . The matter-light interaction creates a periodic potential whose geometry and amplitude are determined by the 
laser configuration and intensity. Both can be controlled in experiments. For instance, using pairs of counter-propagating laser beams 
(Box 1 Fig lb), the lattice potential reads y, z) = V° cos (2kLfT), where is the lattice depth and the laser wavevector. 

In deep lattices, the atoms get trapped at the periodically-arranged minima of the lattice potential (so-called lattice sites). They can jump 
from site to site via quantum tunnelling (with a rate J ), and two atoms interact only in the same site (with an energy U). This physics is 
governed by the Hubbard Hamiltonian, i.e. the discrete version of Hamiltonian (1): 


ft ~ ^ Jtr {o^ jOcrj + h.C.) + ^ V a j Q^ jClo-j + ^ 




3 + a 4 , 

cr,j a-' 


LIq-'jUcj-j 


( 2 ) 






where the sum over (j, l) covers all sites j and their nearest-neighbour sites /, and a a j is the annihilation operator of an atom in site j. 
Hence, ultracold atoms (bosons or fermions) in optical lattices mimic the Hubbard model, which is widely considered in condensed- 
matter physics, for instance to capture the essential physics of electrons in solids. However, in contrast to condensed-matter systems, 
Hamiltonian (2) can be shown to be exact in the limit of deep lattices, low temperature and low interactions 19 . The parameters JVo-j 
and Uo-'O-'j in Box 1 Eq. (2) can be calculated ab initio from the potential Vo-(r) —» Vo-(r) + V^ a (r) and the coupling constant g^ in 
Box 1 Eq. (1) and are thus controllable in experiments. 
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Box 1 Figure 1 | Confining ultracold atoms in magnetic and optic traps, a) Harmonic trapping and laser waveguide (coutesy of 
V. Josse and P. Bouyer). Magnetic coils create a nearly harmonic trapping potential at the bottom of which a degenerate quantum gas, 
surrounded by a cloud of thermal atoms, is formed. A focused laser beam which creates an almost one-dimensional waveguide is also 
represented, b) Optical lattice. The interference pattern of pairs of counter-propagating laser beams form a periodic potential (represented 
here in two dimensions). The atoms are trapped in the lattice sites, but they can tunnel from site to site with a tunnelling rate J and interact 
when placed in the same site with an energy U. 


and numerical results indicate that repulsive interactions decrease the 
localization length before completely destroying localization 77 . Since 
a non-linear term is naturally present in BECs (see last term of Eq. §). 
and can be controlled via Feshbach resonances 45 , transport experi¬ 


ments with interacting condensates are particularly promising to ad¬ 
dress this question. 

A different approach to the interplay of interactions and local¬ 
ization consists in considering a Bose gas at equilibrium in a d- 
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Figure 2 | Effect of interactions in disordered Bose and Fermi gases. 

The gas is described using the single-particle (non-interacting) states !*■,). In 
the presence of disorder, these states, which are localized and distributed 
around in a given volume, are represented by the spheres (in red when they 
are populated). Bose gas: For non-interacting bosons (top, central panel), 
the ground state, \xo), only is populated. Then, attractive interactions (top, 
left panel) tend to contract this state, thus favoring localization. Conversely, 
repulsive interactions (top, right panel) work against localization by populating 
more and more \x v ) states. Fermi gas: In the absence of interactions, a gas of 
N fermions populates the N lowest-energy [y„> states (bottom, central panel). 
Then, each state tends to extend under the action of attractive interactions 
as for maximizing the overlap between different populated |^ v > states (bottom, 
left panel). Conversely, for repulsive interactions, they tend to minimize their 
mutual overlap, then favoring localization (bottom, right panel). 


dimensional box of volume f2 in the presence of interactions and dis¬ 
order (Fig. ^). For vanishing interactions and zero temperature, all 
bosons populate the single-particle ground state, \xo)- Very weak at¬ 
tractive interactions are expected to favor localization by contracting 
the Bose gas, but also induce instabilities for moderate interactions, 
pretty much like for trapped BECs 13 . Conversely, weak to moder¬ 
ate repulsive interactions do not affect much the stability, but work 
against localization, by populating an increasing number of single¬ 
particle states, [yy). Weak interactions populate significantly only the 
lowest-energy states. Since they are strongly bound in rare, low-energy 
modulations of the potential, their mutual overlap is small. The gas 
then forms a Fock state, ]'F) °c Ilv(^v) iVv |0), where b\ is the creation 
operator in state [yy). The population N v of each is determined by the 
competition between single-particle energy e v and interaction within 
each state \x v ). This results in the characteristic equation of state 78 , 
Ng = f M de 2) n (e)(jU - e)P{e), where Da{e) is the density of states and 
Pv = 1 / f dr Ltv(r)| 4 is the participation volume of |^y). This state is 
an insulator with finite compressibility, k = dN/d^, and can thus be 
refered to as a Bose glass 49,50 . It attains particularly interesting fea¬ 
tures in disordered potentials bounded below (i.e. when F(r) > V mm 
everywhere), for which Lifshits has shown 79 that the relevant single¬ 
particle states are determined by large-scale modulations of the po¬ 
tential. Since they are exponentially far apart, the density of state 
is exponentially small, Dn(e) ~ exp[-c(e- Vmin) - ^ 2 ]- As one can 
see, the equation of state is determined by both the density of state 
Dn(e) and the localization via P(e ) in the Lifshits tail, which leads us 
to name this state the Lifshits-Anderson glass 78 . In the opposite limit 
of strong interactions, there are very many populated \Xv), which thus 
overlap, and the above description breaks down. In turn, the gas forms 
an extended, connected (quasi-)BEC of density n(r) = [jj - V(r)]/g, 
which is well described in meanfield approach 8 ". This state is a su¬ 
perfluid. Finally, the intermediate region interpolates between the 
Lifshits-Anderson glass and the BEC regime. Then, the Bose gas sep¬ 
arates in fragmented, forming a compressible insulator (Bose glass). 
The characteristic features of the fragments can then be estimated in 
the meanfield framework 81 . 

The above description is consistent with the idea that even weak 
interactions destroy single-particle localization. In order to gain fur¬ 
ther insight, it is worth noting that in interacting systems, the rele¬ 


vant states are not the single-particle eigenstates, but are of collective 
nature. For interacting BECs, they are Bogolyubov quasi-particles 13 . 
One then finds that, although the ground state is extended, the Bo¬ 
golyubov quasiparticles are localized 47,48,82 . Their localization prop¬ 
erties however differ from those of Schrodinger particles, owing to 
a strong screening of disorder 48 . In ID, the Lyapunov exponent of 
a Bogolyubov quasiparticle reads F(k) = [«S(k)] 2 y(k), where y(k) is 
the single-particle Lyapunov exponent and S{k) = 2(fcf) 2 /(l + 2(fcf) 2 ) 
is the screening function. One thus finds that in the phonon regime 
(k <K 1 /£), the screening is strong and r(k) <K y(k). Conversely, in the 
free-particle regime (k » 1 /£), the screening vanishes and F(k) == y(k). 
Hence, surprisingly, localization can survive in the presence of strong 
mean-field interactions. This poses new challenges to ultracold atoms: 
Not only direct observation of many-body AL, but also possible conse¬ 
quences on quantum coherence, sound-wave propagation or thermal- 
ization process. 

Fermi systems and ‘dirty’ superconductors 

Consider now a Fermi gas, and focus again on the ground state prop¬ 
erties (Fig. ^). In the absence of interactions, the gas of N fermions 
populates the N lowest single-particle states. For low density, short- 
range interactions do not play a significant role as the populated states 
are spatially separated. However, for large-enough density, they do 
overlap. Then, for repulsive interactions, each populated state tends to 
contract to minimize its overlap with other populated states, thus favor¬ 
ing localization. Conversely, for attractive interactions, the populated 
states tend to extend to maximize their overlap, thus favoring delocal¬ 
ization. Hence strikingly, interactions have opposite consequences for 
fermions and bosons. 

Perhaps even more fascinating is the possibility to study ‘dirty’ 
Fermi liquids. Experiments with two-component Fermi gases (e.g. 
6 Li or 40 K), with interactions controlled by Feshbach resonances, 
have already significantly advanced our understanding of the so-called 
BEC-BCS crossover 15,16 . On the attractive side of the resonance and 
for weak interactions, the Fermi superfluid is well described by the 
Bardeen-Schrieffer-Cooper (BCS) theory and formation of spatially 
extended Cooper pairs consisting of two fermions of opposite spins 
and momenta. On the repulsive side, pairs of fermions form bosonic 
molecules, which undergo Bose-Einstein condensation. Although dis¬ 
order should not significantly affect pairing, BCS superfluidity and 
BEC superfluidity are expected to react differently to disorder 83,84 . The 
famous Anderson theorem 85 indicates that disorder should not affect 
very much the BCS superfluid owing to the long-range and overlap¬ 
ping nature of the Cooper pairs. Conversely, disorder should seriously 
affect the molecular BEC, enhancing phase fluctuations. 


Strongly-correlated gases in disordered lattices 

Strong interactions are also very important in various disordered sys¬ 
tems, e.g. superfluids in porous media or ‘dirty’ superconductors. 
Metal-insulator transitions attain a particularly interesting, but not 
fully understood character in lattice systems. In this respect, the Bose- 
Hubbard model, 


H 


- ^ j(d'.di + h.c.) + ^ V) a\cij + 

W) J 




'pl'lM: (4) 


is central in condensed-matter physics 50-52 for it forms a tractable 
model, which captures the elementary physics of strongly interacting 
systems. Hamiltonian (Q) describes bosons, in a lattice with inhomo¬ 
geneous on-site energies Vj , which can tunnel between the sites, with 
rate J , and interact when placed in the same site, with interaction en¬ 
ergy U. Interestingly, this model contains the most fundamental two 
phenomena underlying metal-insulator transitions. They correspond 
to the Anderson transition 5,6 in the absence of interactions ((7 = 0) as 
discussed above, and to the Mott transition 86 in the absence of disor¬ 
der (V = 0). In systems dominated by repulsive interactions, density 
fluctuations, which are energetically costy, are suppressed, and a Mott 
insulator (MI) state, ] v F M i) oc ■ (flj)"|0), is formed. Then, the number 
of bosons per site, n = \ji/U + 1], where [.] represents the integer part, 
is determined and phase coherence between the lattice sites vanishes. 
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Box 2 | Creating controlled disordered potentials. 

In atomic gases, disorder can be created in a controlled way. For instance, the so-called speckle potentials are formed as follows 109 . 
A coherent laser beam is diffracted through a ground-glass plate and focused by a converging lens (Box 2 Fig. la). The ground-glass 
plate transmits the laser light without altering the intensity, but imprinting a random phase profile on the emerging light. Then, the 
complex electric field £(r) on the focal plane results from the coherent superposition of many waves with equally-distributed random 
phases, and is thus a Gaussian random process. In such a light field, atoms with a resonance slightly detuned with respect to the laser 
light experience a disordered potential V(r) which, up to a shift introduced to ensure that the statistical average (V) of l/(r) vanishes, is 
proportional to the light intensity, l/(r) oc ±(|£(r)| 2 - (|£| 2 )), an example of which in shown in Box 2 Fig. lb. Hence, the laws of optics 
allows us to precisely determine all statistical properties of speckle potentials. First, although the electric field £(r) is a complex Gaussian 
random process, the disordered potential F(r) is not Gaussian itself, and its single-point probability distribution is a truncated, exponential 
decaying function, P(V(r)) = e 1 1 Vr| 1 exp(-V , (r)/V , R )0(F(r)/V , R + 1), where -y/( V 2 ) = |V R | is the disorder amplitude and 0 is the 
Heaviside function. Both modulus and sign of V R can be controlled experimentally 31 : The modulus is proportional to the incident laser 
intensity while the sign is determined by the detuning of the laser relative to the atomic resonance (V R is positive for 'blue-detuned' laser 
light 31 ' 38 ' 41,44 , and negative for ‘red-detuned' laser light 37,39 ' 42 ). Second, the two-point correlation function of the disordered potential, 
CSfr) = (V(r)V(0)), is determined by the overall shape of the ground-glass plate but not by the details of its asperities 109 . It is thus also 
controllable experimentally 31 . There is however a fundamental constraint: Since speckle potentials result from interference between light 
waves of wavelength Tl coming from a finite-size aperture of angular width 2a (Box 2 Fig. la) they do not contain Fourier components 
beyond a value 2 k c , where k c = (2 n/A L ) sinter). In other words, C2(2k) = 0 for |k| > k c . 

Speckle potentials can be used directly to investigate the transport of matter-waves in disordered potentials 37-40 . They can also be 
superimposed to deep optical lattices 88 . In the latter case, the physics is described by Box 1 Hamiltonian (2) with V aj a random variable 
whose statistical properties are determined by those of the speckle potential. In particular, V, r j is non-symmetric and correlated from site 
to site. Yet another possibility to create disorder in deep optical lattices is to superimpose a shallow optical lattice with an incommensurate 
period 40 ' 43 ' 45,87 . In this case, V, r j = Acos(2 nj3j + (/>), where A and if> are determined by the amplitude and the phase of the second lattice 
and j8 = fc/fti is the (irrational) ratio of the wavevectors of the two lattices. Although the quantity V, r j is deterministic, it mimics disorder 
in finite-size systems 33 ' 34 ' 89,90 . In contrast to speckle potentials, these bichromatic lattices form a pseudo-random potential, which is 
bounded ( \Vo-j\ < A) and symmetrically distributed. 




Box 2 Figure 1 | Optical speckle potentials, a) Optical configuration, b) Two-dimensional representation of a speckle potential. 


Mis are insulating, incompressible, and gapped as the first excitation 
corresponds to transfer one atom from a given site to another, which 
costs the finite energy U. At the other extreme, when tunneling dom¬ 
inates, the bosons form a superfluid state, | v Psf) x (Z/ o]) |0), with 
normal density fluctuations and perfect coherence between the lattice 
sites. This state is gapless and compressible. 

In the presence of disorder, a glassy phase is formed, which in¬ 
terpolates between Lifshits-Anderson glass for weak interactions, to 
Bose glass for strong interactions 30 . The latter can be represented 
as I'I'bg) 00 [I, (S])” J |0) with nj = [(/* - Vj)/U + 1]. This phase is 
thus insulating but compressible and gapless since the ground state 
is quasi-degenerated, like in many other glassy systems 50-32 . With 
the possibility of realizing experimentally systems exactly described 
by Hubbard models (Boxl), ultracold atoms in optical lattices offer 
also here unprecedented opportunities to investigate this physics in de¬ 
tail, and to directly observe the Bose glass, which has not been pos¬ 
sible in any system so far. Two experimental groups have made the 
first steps in this direction 87 ' 88 . The experiment of ref. 87 applied a 
bichromatic, incommensurate lattice to ID Mott insulators. Increas¬ 
ing disorder, a broadening of Mott resonances was observed, suggest¬ 
ing vanishing of the gap and transition to an insulating state with a 
flat density of excitations. Intensive theoretical studies have been de¬ 


voted to understand these results, using quantum Monte-Carlo 89 and 
Density Matrix Renormalization Group 90 techniques. The results of 
ref. 89 suggest that, in the conditions of ref. 87, one should expect 
a complex phase diagram with competing regions of gapped, incom¬ 
pressible band-insulator, and compressible Bose glass phases. Clearly, 
novel and more precise detection schemes are needed to characterize 
this kind of physics, such as direct measurements of compressibility 51 
or condensate fraction in superfluid, or coexisting superfluid and MI 
phases. The latter has been approached experimentally in ref. 88, 
where disorder-induced suppression of the condensate fraction in a lat¬ 
tice with super-imposed speckle was observed. 

One can also investigate the corresponding Fermi counterparts 
with ultracold atoms. These systems are particularly interesting as 
they would mimic superconductors, even better than bosons. In this 
respect, an outstanding challenge is definitely to understand high-Tc 
superconductors, and possibly important effects of disorder in these 
systems. Consider the two-component (<x e jt.il) Fermi-Hubbard 
Hamiltonian 

H = - Jcr + h-c.) + Vo-j al^o-j (5) 

c r,(j,l) O-J 

+ 'El U 
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For weak interactions, we have a Fermi liquid similar to that discussed 
above. For strong interactions and low temperature, T < U, the Fermi 
gas enters a MI state, pretty much like for bosons, but with a single 
0? = 1) fermion per site (either f or J,). Evidence of vanishing double 
occupancy and incompressibility have been reported recently in Fermi 
Mis 24-25 . Then, in the presence of disorder, various phases could be 
searched for, such as Fermi glasses. At even lower temperatures, spin 
exchange starts to play a role, and a transition from paramagnetic to 
antiferromagnetic insulator phases is predicted for T N ~ 4 J 2 /U in non- 
disordered systems 91 . Interestingly, the interplay of interactions and 
disorder might lead to appearance of novel ‘metallic’ phases between 
the Fermi glass and the MI. Hence, dynamical mean-field theory 92 at 
half-filling predicts that disorder tends to stabilize paramagnetic and 
antiferromagnetic metallic phases for weak interactions. For strong 
interactions however, only the paramagnetic Anderson-Mott insulator 
(for strong disorder) and antiferromagnetic insulator (for weak disor¬ 
der) phases survive. 



Simulating disordered spin systems 

In condensed-matter physics, other important paradigm models where 
disorder induces subtle effects are spin systems, described by the 
Hamiltonian 

«=- H ini § j ■ § i + J i. §,j j ■ § i + J h § ) - § i) - Z h ' ■ S * (6) 

W J 

with either random spin exchange, Jjj, or random magnetic field, h,. 
Ultracold gases can also simulate this class of systems, although not as 
straightforwardly as for Hubbard models. Consider a two-component 
(Bose-Bose or Fermi-Bose) ultracold gas in an optical lattice, as de¬ 
scribed by Box 1 Hamiltonian (2). In the strongly-correlated regime, 
the couplings between the particles can be understood as exchange- 
mediated interactions between composite (bosonic or fermionic) par¬ 
ticles. One can then map Box 1 Hamiltonian (2) onto Hamiltonian (^) 
with fictitious spins encoded in combinations of the annihilation and 
creation operators of the composite particles: S ^ = (Aj + At)/2, 
Sj' y = (Aj - At)/2i, and S^ = 1/2 - A^A, which indeed have com¬ 
mutation relations of spins 2 , and complicated but analytical functional 
dependence of Jjj and h ; on the parameters of Box 1 Hamiltonian (2). 
In the presence of disorder, these parameters are random 54 59 and one 
can reach various limiting cases corresponding to Fermi glass, quan¬ 
tum spin glass and quantum percolation 58 . 

Particularly promising is the possibility of simulating spin 
glasses 58 (Fig. HI), for which only the exchange term, Jjj, is randomly 
distributed. The phase diagram of (even classical) spin glasses, which 
is not known yet, is an outstanding challenge in condensed-matter 
physics. The nature of spin glasses is still debated and there exist 
competing theories: The Parisi replica symmetry breaking 56 and the 
Nelson-Huse droplet model 57 . Ultracold atoms might contribute to the 
resolution of this issue, not only on the classical level but also on the 
quantum level since they offer original ways of performing quenched 
averages. Importantly with a view towards testing the replica theory, 
overlap between two spin configurations between two (or more) repli¬ 
cas can be measured directly by preparing a pair of 2D lattices with the 
same realization of disorder 93 . Quenched averages for systems with bi¬ 
nary disorder can also be simulated by replacing the classical disorder 
variables by quantum 1/2-spins, and preparing them in a superposition 
state 94 . 

Yet another fascinating possibility is to simulate various random 
field-induced order (RFIO) phenomena in systems with continuous 
symmetry, such as BECs or XY-spin models with 1/(1) symmetry, or 
Heisenberg models with S (7(2) symmetry 54 55 . A prototype model 95 , 
is the 2D-XY version of Hamiltonian (n) with fixed exchange Jjj 
but random field h^. In the absence of disorder, symmetry leads to 
strong fluctuations, which suppress long-range order, according to the 
Mermin-Wagner-Hohenberg theorem. Disorder distributed in a sym¬ 
metric way suppresses ordering even more. Surprisingly however, dis¬ 
order that breaks symmetry might actually favor ordering. This model 
can be implemented within Bose-Bose mixtures 54,55 where random 
uniaxial h ; can be implemented using two internal states of the same 


Figure 3 | The spin glass problem. An assembly of spins located at the 
nodes of a cubic lattice interacts according to Hamiltonian (|| where the ex¬ 
change term Jjj only is randomly distributed, and can be either ferromag¬ 
netic (blue bonds) or anti-ferromagnetic (red bonds). The ground state of the 
system corresponds to the spin configuration that minimizes the total energy. 
The inherent complexity of spin glasses results from frustration which appears 
when the topography of ferromagnetic and anti-ferromagnetic bonds make im¬ 
possible to fulfil the local constraints all together. In some plaquetes of four 
sites, local minimization is easy, for instance when all bonds are ferromagnetic 
(left disk) or anti-ferromagnetic (central disk). In some others, it leads to frus¬ 
tration, for instance for odd numbers of ferromagnetic and anti-ferromagnetic 
bonds (right disk). In the latter case, at least one spin is frustrated, that is its 
spin orientation is not unique. Hence, frustration is at the origin of a manifold of 
metastable states which corresponds to configurations with similar energies. 

atom, coupled via a random Raman field, (zf2(r)T / 1 (r) tv P2(r) + h.c. 
In order to break the continuous symmetry, one uses a Raman cou¬ 
pling with constant phase, but random strength. In lattice systems, 
RFIO shows up but is limited by finite-size effects, even in very large 
systems 54 . In this respect, ultracold atoms offer an alternative and fruit¬ 
ful route. Indeed, RFIO turns out to be particularly efficient in two (or 
multi-) component BECs in meanfield regime, where the energy func¬ 
tional reads A E ^ dr ;?[(S 2 /4m)(V0) 2 +/iO(r) cos(0)], with n the atomic 
density and Q{ r) = £fi(r) - 9j(r), the phase difference between the two 
BECs. This is the continuous counterpart of the 2D-XY model. Then, 
RFIO manifests itself as a fixed 6( r) = ± zr/2, and thus allows to con¬ 
trol the relative phase between the components 55 . This is a striking 
example where ultracold atoms can be used not only to simulate clas¬ 
sic models, but also offer new and fruitful viewpoints to fundamental 
issues. 

Further directions 

As the reader has probably noticed, we both are very enthusiastic 
about the future development of the field of disordered quantum gases, 
and probably would like that any interesting direction can be pursued. 
Limited size of the present review has not allowed us to discuss them 
all, but let us briefly mention another. 

Two-component (Bose-Bose or Bose-Fermi) mixtures offer an al¬ 
ternative method to create disorder in optical lattices, namely by 
quenching one component in random sites, so as to form a background 
of randomly-distributed impurities 94,96 . Theoretical analysis using 
Gutzwiller method confirms the appearance of incompressible MI and 
partially compressible Bose glass phases 97 . The idea of freezing the 
motion of the second species to form random impurities (i.e. classical 
disorder) can be generalized to freezing of any quantum state 98 . In this 
case the system does not involve any classical disorder, but neverthe¬ 
less localization occurs owing to quantum fluctuations in the frozen 
state of the second species. 

One can even relax the freezing condition and consider say two 
bosonic species, one of which tunnels much slower than the other, 
forming a quasi-static disorder. In a large region of parameters (for re¬ 
pulsive inter-species forces), the ground state corresponds to full phase 


6 































segregation. In practice it is marked by a large number of metastable 
states in which microscopic phase separation occurs, reminiscent of 
emulsions in immiscible fluids". Such quantum emulsions are pre¬ 
dicted to have very similar properties to the Bose glass phase, i.e. 
compressibility and absence of superfluidity. Such quasi-static or even 
time-dependent disorder effects have been suggested to underlie the 
quite large shift of the SF-MI transition in Bose-Fermi 100-101 and Bose- 
Bose 102 mixtures. This issue was quite controversial and the most re¬ 
cent work suggests that, while indeed the fermions tend to localize the 
bosons for attractive boson-fermion interactions, higher Bloch bands 
play a significant role 10 ’ -105 . 
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